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^ Abstract 

We consider fermions in a zero-temperature superconducting anti-de Sitter domain wall 
solution and find continuous bands of normal modes. These bands can be either partially 
filled or totally empty and gapped. We present a semi-classical argument which approxi- 
mately captures the main features of the normal mode spectrum. 



November 2009 



Contents 



1 Introduction 



1 



2 The bosonic background 



3 



3 A semi-classical argument 



4 



4 Bands of fermion normal modes in AdS± 



9 



5 Discussion 



1 Introduction 

The Abelian Higgs model in anti-de Sitter space [U [2] is the most straightforward way to 
realize superconducting black holes in string theory [31 HI [5] . At the level of supergravity 
calculations, it is a purely bosonic model: the lagrangian involves gravity, an abelian gauge 
field, and a complex scalar which is the order parameter for superconductivity. Already in 
PP it was suggested that the complex scalar should be the dual of an operator which destroys 
a Cooper pair. However, the dual field theory typically involves both fermions and bosons. 
For example, in the construction of [3], the dual of the complex scalar is a sum of a fermion 
bilinear and a scalar trilinear. One might imagine a case where the field theory has no scalars 
(or, at least, none charged under the global symmetry dual to the U(l) gauge symmetry of 
the bulk that gets spontaneously broken). Although no such example has been exhibited 
explicitly, there also isn't any argument we know of that there can't be one. 

With fermions in the field theory, one can certainly consider color singlet operators dual 
to fermions in the bulk: for example, operators schematically of the form trX n A with n > 1 
could be dual to spin- 1/2 fermions, and tr X n \XX with n > could be dual to either spin- 1/2 
or spin-3/2 fermions, depending on how indices are contracted. (More properly, one must 
anticipate that bulk fermions are dual to a linear combination of operators with an odd 
number of field theory fermions.) It's difficult to have a bulk fermion dual to an operator 
with one field theory fermion and nothing else, because it's individual fermions in the field 
theory aren't gauge singlets. Calculations with fermions [6J [Zl El E] have focused on an 
properties of the two-point function of fermionic operators in the background of a Reissner- 
Nordstrom-anti-de Sitter black hole (hereafter RNAdS), in particular a singularity at finite 
momentum which has been argued to be evidence of non-Fermi liquid behavior. 
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In this paper we wish to study the behavior of spin- 1/2 bulk fermions in response to the 
superconducting AdS± domain wall solution of [5] , a zero-temperature geometry whose finite- 
and low-temperature limits were previously studied in [I]. This domain wall solution has 
two advantages over the four-dimensional RNAdS solution: 1) RNAdS is unstable toward 
superconducting instabilities, and the domain wall solution is plausibly the endpoint of the 
evolution of one such instability; 2) RNAdS has macroscopic entropy at zero temperature 
which has not found a satisfactory explanation, but the domain wall solution has no entropy 
at all, at least in the classical supergravity approximation. A technical disadvantage of 
the domain wall solution is that it is known only numerically, so to study fermions we will 
have to solve differential equations whose coefficients are known only numerically. Based 
on studies including [101 [HJ [121 H3], we expect that domain wall solutions with at least 
Lorentz invariance in the infrared are fairly generic ground states of the Abelian Higgs model, 
although Lifshitz scaling in the infrared is another possibility. Domain walls with conformal 
invariance in both the ultraviolet (UV) and infrared (IR) are under the best theoretical 
control, since curvatures can be made small everywhere. But we anticipate that arguments 
presented here could be extended to more general domain wall solutions, and the conclusions 
might be fairly similar when there is emergent Lorentz symmetry in the infrared. 

Ideally, the fermions we should study in the AdS± domain wall should be the ones present 
in maximal d — 4 supergravity (or some consistent truncation thereof). There are two reasons 
we do not do this. The first is simplicity: The quadratic actions for fermions in these theories 
are complicated, and there appears to be significant mixing between the gravitini and spin- 
1/2 fermions. The second is flexibility: We will consider different values of the charge and 
mass of the fermion, and we will see that the final results significantly depend on the choice 
of these parameters. 

The remainder of this paper is structured as follows. In section [2] we briefly review the 
domain wall background. In section [3] we present a semi-classical argument which focuses 
attention on a compact region of phase space and suggests that the normal modes will lie 
approximately along segments of hyperbolas. In section [4] we numerically solve the Dirac 
equation for a charged fermion in the AdS± domain wall solution, finding one or more bands 
of fermion normal modes. We conclude with a discussion in section [5l 

When this paper was nearing completion, we received [H], which studies fermions in a 
related geometry. We were also informed by D. Vegh that he, T. Faulkner, G. Horowitz, 
J. McGreevy, and M. Roberts are also working on related topics. 
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2 The bosonic background 



The bosonic lagrangian that is the basis for the domain wall background that we will study 
takes the form 



jC — R — —Fj,,, — — 
4 2 



(^77) 2 + sinh%^-^^ 



+ — cosh 2 ^(7 — cosh 7]) 
L 1 2 

where rj is a real scalar, 9 is a real pseudoscalar, and we use mostly plus metric conventions. 
We refer the reader to the literature [5j [10] for details on the construction of the domain wall 
solution, but summarize some of its properties for convenience. The domain wall geometry 
takes the form 

2B(r) J 2 

L J h[r) 

and has non-zero i](r) and gauge field A^dx^ 1 = Q(r)dt. We can always choose coordinates 
where B = 0, and will make this choice in numerical computations but will keep B general in 
formulae. The scalar potential in has two extrema, at rj = and 771R = log(3 + 2 3 / 2 ), and 
associated to these are AdS^ solutions with radii of curvature Luv — L and Ljr = \/3L/2, 
respectively. The domain wall interpolates between these two AdS geometries. We can 
choose coordinates such that as r —>■ — 00 

A ~ — h ~ 1 r] ~ r] lR $ ~ , (3) 

Li — 



IK 



while as r — » +00, 

r 

v4 ~ /l ~ h\jy 7] ~ $ ~ $uv , (4) 

where fouv and $uv can only be determined numerically and are given by /iuv ~ 14.249 
and $uv = 9.328. The solution is shown in Fig. [TJ At a constant r slice the metric is, 
up to a constant rescaling, the Minkowski metric with an effective "speed of light" given 



by v (r) = yjh(r). The value of v(r) can be changed by a redefinition of the t and x 
coordinates, but ratios of v at different values of r are invariant. In our choice of coordinates, 
f IR = v(— 00) = 1 and Vyy = v(+oo) = 3.775. 
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Figure 1: The metric, gauge field and scalar field for the domain wall solution in M-theory. 

3 A semi-classical argument 

The fermionic lagrangian that we will consider is 

£ f = ^(r^-m)^, (5) 

where D^ip = (d^ + \u^T pa — iqAA ip and is the spin connection. We employ the 
same conventions on fermions in AdS^ as in [6], except that we use /i for a curved space 
index and fi for a tangent space index. We always assume m > 0. 

Much can be learned from the asymptotics of solutions to the Dirac equation following 
from ([5]) in the UV and IR asymptotic regions. These asymptotics are simple because all 
one is describing is a free fermion propagating in empty anti-de Sitter space with a constant 
electrostatic potential $. Our discussion here will be somewhat heuristic, relying on replacing 
the Dirac equation by an on-shell relation which is only justified in a geometric optics limit. 
We will give a more precise treatment in section |4j 

The on-shell constraint implied by the Dirac equation is 

gr(k» - qA,)(K - qA v ) + m c 2 ff = . (6) 

Here A^ = ($, 0, 0, 0) and k p = (—u, £4, k 2 , k r ). Without loss of generality one can use the 
rotational symmetry in the x 1 -x 2 directions to set k 2 = and k\ = k > 0. The frequency uj 
and the momentum k are definite real numbers when we choose the wave-function to take 
the form e - tuJt + lkxl times some function of r. The radial momentum varies as a function of r. 
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The effective mass is the bare mass m? plus some contributions from the spin connection 
and curvatures. In components, ^ becomes 

+ + k 2 + e 2 ^k 2 + e 2A m 2 eS = . (7) 

We next observe that the last term is suppressed in the infrared compared to that first two, 
unless m 2 s increases quickly in the infrared. This is impossible if the infrared geometry is 
anti-de Sitter space: indeed, in that case, m 2 s is constant in the infrared. In geometries where 
only Lorentz symmetry is recovered in the infrared, it is possible that m 2 s does increase very 
quickly, either because of the contributions from the spin connection and curvatures, or from 
a dependence of m on some scalar field which diverges as one passes to the extreme infrared. 
Barring such a circumstance, we see that the sign of k 2 is the same in the infrared as the 
sign of (u + g$) 2 //iiR — k 2 . Limiting ourselves now to the case where $ — > 0, we see that k 2 
is non-positive when 

ImI / — 

< k where v ir = y hiR . (8) 



M 

VlR 

Now, if k r were real an non-zero, then its value would be essentially the radial momentum of 
the fermion. As an extension of standard boundary conditions at a black hole horizon, one 
may reasonably require this momentum to be infalling (that is, toward the infrared). But it 
makes no sense for the fermion to be falling down ever further into the infrared when its state 
is a normal mode. Instead, k r should be either or imaginary, so that the wave-function 
of the fermion can decay as one passes further into the infrared. If the infrared geometry is 
anti-de Sitter, then the decay is very rapid: exponential in |fc r |e~ A , or power-like in e~ A if 
k r = 0. In a more general setting, the right requirement is for the fermion wave-function to 
be normalizable in the infrared, which presumably rules out oscillatory and singular behavior 
while allowing a more regular solution similar to the exponentially damped solution in anti- 
de Sitter space. The upshot of this discussion is that ^ is one requirement that a normal 
mode must satisfy. 

The requirement that u = for a fermionic normal mode in RNAdS comes from reasoning 
rather similar to what we just went through: only for this value of u can one avoid oscillatory 
behavior in the infrared indicating that the fermion is falling into the black hole. In fact, 
since h — *■ at the horizon of the RNAdS geometry, one has i>xr = 0, so one can still use 
w — vm\k\ as the condition that restricts the possible values of u and k based on the infrared 
dynamics. 

In the extreme ultraviolet, the last term of ^ dominates over the first two, so one can 
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conclude that k 2 has the sign opposite m e 2 ff , namely negative. So oscillatory behavior is 
impossible. Though our reasoning here is non-rigorous, the conclusion that solutions are 
non-oscillatory is correct: all components of the fermion wave-function must have power-law 
behavior in e~ A , no matter what to and k are. Now, it is hard to see how one would find 
a normal mode if k 2 were always negative, because then either the wave-function of the 
fermion would tend to increase monotonically toward the ultraviolet or else monotonically 
toward the infrared, neither of which would be consistent with normalizable behavior. So, 
as a necessary condition beyond d8|), it is reasonable to expect that 



sup k 2 = sup { ( U + «^ 2 - k 2 - e 2A m 2 cS \>0. (9) 



r 



In principle, knowing a bosonic background, one can use (|9| to obtain an upper bound on 
values of k where normal modes can exist, as a function of uo. In practice this is laborious. 
We prefer to use instead the related condition: 



fuv 



uv 



>k, (10) 



where t>uv = V^Hjv- (|9j) and (10) are equivalent in the limit where cu/fuv and k are large 
compared to m e $. So (10), which is simple to apply, is really the large frequency, large 
wave- number limit of the better justified condition (|9]). 

To summarize: A heuristic and inexact treatment of asymptotics in the UV and IR 
regions of a domain wall geometry leads to the expectations that normal modes will be 
confined to what we will term the "preferred wedge," namely 

M< fc< + (ii) 

VlR V V y 

In other words, the projection of the gauge-invariant momentum k^ — qA^ onto boundary di- 
rections should be timelike in the ultraviolet, in order to produce non-monotonic dependence 
on r, and spacelike in the infrared, in order to ensure that the fermions aren't falling into 



the far infrared. The first inequality in (11) (from the infrared constraint) can be expected 
to be more reliable than the second inequality (from the ultraviolet behavior). The preferred 
wedge is compact, as figure [2] shows. 

Having established our expectations for the region of the u-k plane where normal modes 
may appear, the next obvious question is what the pattern of normal modes will be in this 
region. First, if there are normal modes at all, it is reasonable to think that they form 
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Figure 2: The green wedge is where (11) holds. This is the region where our heuristic ge- 



ometric optics arguments indicate that one might find fermion normal modes. The grey 



and black curves are an approximate depiction of the hyperbolas (13), representing an ap- 



proximate WKB treatment of where normal modes lie. Only the black parts of the curve 
correspond to actual normal modes; the grey parts are where normal modes might have been 
if the green wedge had been larger. 



one-dimensional families. The reason is that individual components of the Dirac spinor 
satisfy a second-order differential equation, so there are two linearly independent solutions. 
Demanding that the one that is normalizable in the infrared is also normalizable in the 
ultraviolet amounts to a single real condition on the two real parameters u and k. So indeed 
one expects one-dimensional families. 

A Bohr-Sommerfeld estimate of the location of normal modes takes the schematic form 



dr k r = n(n + v) 



(12) 



7lH.(w,fc) 



Here ruv and are the radii where k r = is a solution to 0; k r is the positive root of 
^ for r IR < r < ruv; n is a non-positive integer; and a standard expectation for v (the 
Keller-Maslov index) is 1/2. Passing to a further approximation, one can consider the case 



where the integral ( 12 ) is dominated by radii where the geometry is approximately AdS- 



uv 
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and the gauge potential is approximately $uv- Then the integral (12) should be roughly 



proportional to \l k 2 , with a constant of proportionality that depends on details 



of the bosonic background and so cannot readily be computed. Using (12), we are led to an 
approximate condition 



k = (Y.n + Y 2 ) 2 (13) 



for the position of the n th band of normal modes. Y\ and Y 2 are constants and n is again a 
non-positive integer, which is the number of nodes in the fermion wave-function. (Different 



components of the Dirac spinor could have bands of normal modes described by (13) with 
different Y 2 but the same Y\. Associating n with the number of nodes can only be expected 
to work when one is focusing on the normal modes involving a particular component.) Evi- 
dently, we should expect only finitely many bands of normal modes, because the hyperbola 



(13) doesn't intersect the region |u;| < v\^k when n is large. 



Since the hyperbolas (13) do not intersect the UV boundary \u — q$\jv\ = v\jyk of the 
preferred wedge region for normal modes, the topological prediction we get out is that the 
bands begin and end on the IR boundary, \u\ = Vj&k. It is possible for a band to both begin 
and end on the upper branch of the IR boundary, i.e. the one with u > 0. Herein lies an 
interesting possibility: for such a band, there would be a minimum energy required to add a 
fermion belonging to it (unless the band dipped below the u = axis and then rose back up). 
This is obviously an attractive feature if one is to make contact between these holographic 
models of superconductors and real-world superconductors. However, we are obliged to add 
that it is equally possible for bands of normal modes to cross the u = axis one or more 
times. Then the natural configuration to consider is the one where the u < normal modes 
are populated, while the u > modes aren't. This results in a Fermi surface in the bulk, 
and apparently a gapless state for the fermions. Two considerations might change this state 
of affairs: 1) The fermions interact gravitationally and so presumably have an attractive 
channel that causes them to superconduct in the bulk through the usual mechanism of 
forming Cooper pairs; and 2) The gas of fermions back- reacts on the bulk, as in [To] . 

In the next section, we will exhibit a case where most bands cross the u> = surface at 
least once, and another case where there is a single band that does not cross the u — axis 
and therefore has a gap in the sense explained in the previous paragraph. It is interesting 
indeed to inquire what behavior the fermions of maximal gauged supergravity exhibit in the 
domain wall geometry summarized in section [2j gapped or ungapped? 

As a final aside, it is easy to imagine cases where there are no normal modes: the con- 
dition (m) could be impossible to satisfy, or it could be that Y 2 in (13) is too large for any 



8 



of the would-be normal modes to intersect with permitted infrared boundary conditions. 
This doesn't seem quite to be a situation meriting the term "gapped behavior," because 
there isn't a definite energy that one can add to the system to obtain a long-lived fermionic 
quasi-particle. But perhaps it is reasonable to expect that in a strongly interacting supercon- 
ductor at zero temperature, unpaired fermions can never be stable. Normal modes could be 
altogether absent in ground states with emergent conformal or Lorentz symmetry, where the 
preferred wedge is (modulo caveats already discussed) really a triangular wedge as shown in 
figure [2] or they could be altogether absent in ground states with Lifshitz symmetry, where 
t>iR = and the preferred wedge shrinks to a line segment along the uo = axis; or, indeed, 
they could be altogether absent in an RNAdS vacuum: this happens if q is too small. It 
would be interesting to study the spectral measure of the fermion two-point function in sit- 
uations where there are no normal modes. It could be that a ridge similar to the one found, 
for example, in [BJ, would remain, suggestive of unstable fermionic quasi-particles. 



4 Bands of fermion normal modes in AdS± 

Let us now describe the numerical computation of the normal modes in more detail. It is 
convenient to choose a basis of V matrices such that V- is diagonal, such as 



1 
-1 



y 
V 



z = 0,l,2, 



(14) 



where 7* are gamma matrices in three dimensions given by 7 = i<72, 7 1 = 01, and 7 2 = 03. 
As argued above, symmetry allows us to take the spinor to be of the following form: 



ip{t,x,r) = e- iuJt+tkxl u(r) = e ~^ ikxl 



fuf(r)\ 
ut(r) 
u i( r ) 

\ u 2 ( r V 



(15) 



The Dirac equation can then be written as 



6hA' + h' r 
— — 1 - — m 

4Vh 



u = 



(16) 



Let us now consider the asymptotic behavior of solutions to (16). We are interested in 



regular and purely infalling solutions, which implies that in the infrared (corresponding to 
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u±(r) « t^e- 2A A:i ±mLlH (« m e-^ « m = L my Jk*-u*. (17) 
Here a = 1, 2, is a modified Bessel function, and the £7"^ are constants that satisfy 



U, = _U^ = h+u; 
U+ Ui V fa-u 

and are otherwise arbitrary. In what follows, we will set = C/+ = 1. 

In the ultraviolet, corresponding to r — > +oo, the most general solution is of the form 

«J = e"^ i_ mLuv («:uve- A ) + D^ A I^ +mLuv [^e~ A ) , (19) 

where Kyy = k\— M^guv) anc j j g a mo dified Bessel function. As before, (16) allows you 



to solve for C a , D a algebraically in terms of , D£ , but the exact form of this relation will 
not be important. 



Expanding Ql 9p at large r, we see that 

(20) 



Ua a I 2 J r(i-mL uv ) a I 2 J rd-rnLw) ' 



According to the standard AdS/CFT prescription, we identify with the source and 
with the response of a boundary fermionic operator with conformal dimension A = 
| + mLuvQ 

Normal modes are then solutions of ( 16 ) that in the infrared satisfy ( 17) and have = 0. 
They correspond to poles of the retarded Green's function. To make this connection more 
explicit, consider the full fermionic action 

S f = i J d 5 xij(T»D^ - m)V> + S hdy (21) 

where Sbd y) which we neglected up to now, is a boundary term that does not contribute to 
the equations of motion. It is necessary to have a well posed variational problem and gives 



1 For < m < 1/2, it is also legitimate to identify D+ with the source and with the response, but we 
will in this paper stick with the identification valid for all m. 
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the only nonzero contribution to the onshell action. We choose 

S'bdy = -i J d 4 x V-gg rr , if>± = ^ (1 ± T-) ip , ^ 

r=l/e 

where e is a positive quantity to be taken to zero after functional derivatives are taken. By 



taking the appropriate functional derivatives of (21 ), we obtain the retarded Green's function 

m\ 



//t uv \ mL uv r (I - mLuv) 



;5 o 



G -=m " rtfr^yi Cl (23) 



and we see that zeros of C£ are poles of Gn while zeros of C± are poles of G22- From (16) 



and the form of the T matrices (14) it is clear that the four components of the spinor u 
couple to each other only in pairs: and u\~ mix, and so do u\ and u\. This means we 
are free to consider only u\ and nonzero and look for zeros of or consider only u\ 
and Ui nonzero and look for zeros of G% . 

Before discussing our numerical results let us comment on the meaning of the "preferred 
wedge" (11). It is clear that (J8]) translates^] to «j R > which implies that u is not oscillatory 



in the IR. Condition (10) on the other hand, translates to «uv < which implies that u is 
oscillatory in the UV. This is in agreement with the geometrical optics arguments that led 
us to formulate these conditions in section El 

We now have the necessary ingredients to find normal modes numerically: we can solve 



(16) with initial conditions given by (17) for some large negative r. We can then determine 



the coefficients by fitting the numerical results to (19) for large positive r and vary the 
parameters until we find zeros of the C+. As expected, we find a continuous set of normal 
modes inside the "preferred wedge". As you increase q, the "preferred wedge" grows and 
accommodates more bands. 

As a specific example consider m = and qL\jy = 10, for which the bands are shown 
in Fig. [3] We find an abundance of normal modes, with poles of G22 alternating with poles 
of Gn. The number of nodes in the wave functions (some of which are shown in Fig. [4]) 
is greater for bands that are closer to the origin. The bands are well approximated by 



hyperbola, although (13) does not seem to capture their shape very well. Generalizing (13) 



we can fit the bands to 

(uj + q$) 2 



, ~ = ml s , (24) 

v 



2 Rcmember we are using units where «ir = 1. 



11 



qL uv =10, m=0 



> 
k3 



is 



i 



■ / / / 



/ 



kLu\ 



Figure 3: Fermion normal modes in the AdS± domain wall for m = and qLyy = 10. The 
black lines mark the boundary of the "preferred wedge" (11). The red lines correspond to 



normal modes where and are nonzero (poles in G22) while the blue lines correspond 
to normal modes where the u\ and u{ are nonzero (poles in Gn). The gray dot-dashed lines 
are one parameter fits and the black dashed lines are three parameter fits to hyperbola. The 
red and blue dots mark the location of the normal modes shown in Fig. |4j 



where we treat m e g as an arbitrary fitting parameter. We can either use the extreme UV 
values of $ and v and do one-parameter fits (shown as gray dot-dashed lines in Fig |3| or 
treat $ and v as arbitrary fitting parameters also and do three-parameter fits (shown as 
black dashed lines in Fig pi). In the former case, we find that m e s increases as a small power 
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of the band number. In the latter case (unsurprisingly) we obtain better fits, and we note 
that the best fit values of $ and v approach the expected UV values when you consider 
bands farther from the origin. 




Figure 4: The wave functions for two fermion normal modes of the AdSi domain wall for 
m = and qL\jy = 10. We note that in the Km. > region, we can always choose the w+ to 
be purely real, in which case it follows that the u~ are purely imaginary. The plot on the left 
shows the real part of u\ (blue) and the imaginary part of (green) corresponding to the 
blue dot in Fig. [3j The plot on the right shows the real part of u\ (red) and the imaginary 
part of u{ (cyan) corresponding to the red dot in Fig. pi 



In Fig. |3j all the bands cross the u = line and are therefore ungapped. As you decrease 
q, the "preferred wedge" shrinks and fewer bands are present, but for m = at least one of 
them always seems to intersect the to = line. We can argue that this is always the case if 
we note that for small enough u and k, the Green's function should approach its AdS value 

m 

G R = f (mi IR ) ^--S ■ k oc ° ) , (25) 

V vKy 

where k 2 = k 2 — u 2 and f{mLm) is some unimportant constant. Thus, for m — 0, there is 
always a divergence at the origin that signals a normal mode with (u, k) = and at least 
one of the bands is ungapped. For positive m however, this divergence is gone and we expect 
the possibility of gapped bands. In fact, for gLuv = 3/2 and mL ]R = 1, we find a single 
gapped band (see Fig. [5]). 

Finally, we note that it is also possible to choose the parameters in such a way that there 
are no normal modes at all. If qLjjy is small enough, this seems to be the generic behavior. 
For instance, taking m = and qL\jy = 1/8, we find no normal modes. 
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Figure 5: Fermion normal modes of the AdS± domain wall for mLin = 1 and qLjjy = 3/2. 



The black lines mark the boundary of the "preferred wedge" (11). The blue line corresponds 
to normal modes where the u\ and u{ are nonzero (poles in G n ). Notice it never intersects 
the uj = line (red). So this is a gapped band. 

5 Discussion 

We studied charged Dirac fermions in an AdS^ domain wall solution and found that generi- 
cally they exhibit continuous bands of normal modes. This is in stark contrast to fermions in 
RNAdS, which exhibit one or more isolated normal modes at uj = and finite k [HI [7J El E] . 
However, we note that the same semiclassical argument that led us to expect continuous 
bands in the domain wall geometry can be used to argue for isolated normal modes in 
RNAdS. In fact, since RNAdS has an horizon at which the metric function h vanishes, it can 
in some sense be thought of as degenerate domain wall for which vi& = 0. The condition on 
the location of the normal modes given by([8| then degenerates to uj = and, choosing units 



where t>uv — 1> the remaining condition (10) becomes |<7$uv| > k. The "preferred wedge' : 



thus becomes a "preferred segment" and the same parameter counting that led us to expect 
continuous bands for the domain wall geometry tells us that we should find a discrete set of 
normal modes at uj = 0, which is precisely what happens for RNAdS. 

Since in our conventions, energy equal to the Fermi energy is equivalent to frequency 
equal to zero, the natural zero temperature configuration is one where the normal modes 
with uj < are occupied and those with uj > are unoccupied. Therefore, a particularly 
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important feature of the bands is whether they cross the u = axis. We found this feature 
depends on the choice of the charge and mass of the fermion. For zero mass fermions, the 
bands seem to always cross the u = axis and the resulting configuration would seem to 
have a Fermi surface. However, to better understand this Fermi surface it would be necessary 
to consider the back-reaction of the fermions on the bulk geometry, which might conceivably 
destroy the Fermi surface. For positive mass fermions, with a suitable choice of the charge, 
we found an example of a single band that does not touch the u — axis, i.e., the band is 
gapped. Such a feature is desirable in making comparisons with real-world superconductors 
at zero temperature, and it may be fairly generally achievable when there is at least emergent 
Lorentz symmetry in the infrared. 

Because our bosonic background is embeddable in M-theory jH |5], it would be interest- 
ing to redo our calculations using the quadratic fermion action of maximal gauged d — 4 
supergravity. It appears to be non-trivial to diagonalize this quadratic action, but the ad- 
vantage is that one should in principle be able to understand the operators dual to the 
fermions. It would also be interesting to redo our calculations in AdS$, where we expect 
qualitatively similar results, simply because the semi-classical arguments of section [3] don't 
depend on dimension. Preliminary numerical studies yielded results in agreement with these 
expectations. 
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